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We present simultaneous X-ray, far-ultraviolet, and near- ultraviolet spectra 
of the Seyfert 1 galaxy NGC 7469 obtained with the Chandra X-Ray Observa- 
CN ! tory, the Far Ultraviolet Spectroscopic Explorer, and the Space Telescope Imaging 

\ Spectrograph on the Hubble Space Telescope. Previous non-simultaneous obser- 

vations of this galaxy found two distinct UV absorption components, at —560 and 
■^j- ! —1900 km s _1 , with the former as the likely counterpart of the X-ray absorber. 
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We confirm these two absorption components in our new UV observations, in 
which we detect prominent O VI, Lyct, N V, and C IV absorption. In our Chan- 
dra spectrum we detect O VIII emission, but no significant O VIII or O VII 
absorption. We also detect a prominent Fe Ka emission line in the Chandra 
spectrum, as well as absorption due to hydrogen-like and helium-like neon, mag- 
nesium, and silicon at velocities consistent with the —560 km s _1 UV absorber. 
The FUSE and STIS data reveal that the H I and C IV column densities in 
this UV- and X-ray- absorbing component have increased over time, as the UV 
continuum flux decreased. We use measured H I, N V, C IV, and O VI column 
densities to model the photoionization state of both absorbers self-consistently. 
We confirm the general physical picture of the outflow in which the low velocity 
component is a highly ionized, high density absorber with a total column den- 
sity of 10 20 cm -2 , located near the broad emission line region, although due to 
measurable columns of N V and C IV, we assign it a somewhat smaller ionization 
parameter than found previously, U ~ 1. The high velocity UV component is of 
lower density, logiV = 18.6, and likely resides farther from the central engine as 
we find its ionization parameter to be U — 0.08. 

Subject headings: galaxies: active — galaxies: individual (NGC 7469) — galaxies: 
Seyfert — quasars: absorption lines — ultraviolet: galaxies — X-ray: galaxies 

1. Introduction 

Absorption edges are visible in the X-ray spectra of about one half of all low redshift 
AGN (Reynolds 1997; George et al. 1998; Crenshaw, Kraemer, & George 2003). Ensemble 
studies of the warm absorber population in Seyfert galaxies (Blustin et al. 2005) indicate, 
from estimates of the radial positions of the absorbers, that these systems are likely predom- 
inately outflows from the dusty torus surrounding the accretion disk in each object, while 
those seen in PG quasars are more likely to arise in an accretion disk wind. All of these ob- 
jects also show high-ionization absorption lines in their UV spectra (Crenshaw et al. 1999). 
This suggests a connection between the two phenomena, if not a complete identification of 
one with the other. The position of the absorbing complexes in velocity space, generally 
blueshifted with respect to the systemic redshift of the parent AGNs, the variability of the 
absorption (Kriss et al. 1995; Shields & Hamann 1997; Crenshaw & Kraemer 1999; Cren- 
shaw et al. 2000; Kraemer et al. 2001, 2002; Gabel et al. 2003a) and the non-unity continuum 
and/or broad emission line region (BELR) covering fractions determined for many UV ab- 
sorbers (Kriss et al. 2000b; Gabel et al. 2003b; Kraemer, Crenshaw, & Gabel 2001; Kraemer 
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et al. 2001; Kraemer et al. 2003; Brotherton et al. 2002; Crenshaw et al. 2003) point to an 
origin for the gas that is intrinsic to the AGN. AGN outflows are a potentially significant 
source of feedback energy (Granato et al. 2004; Scannapieco & Oh 2004), influencing the 
luminosity-temperature relation for the intracluster medium of galaxy clusters by steepening 
it on the scale of galaxy groups (Cavaliere et al. 2002), as well as providing metals to and 
governing the formation of galaxies and later generations of AGNs from the intergalactic 
medium (Adelberger et al. 2003). 

Intrinsic absorption in many AGNs has been studied using high spectral resolution UV 
and X-ray observations: Mrk 279 (Scott et al. 2004; Gabel et al. 2005); Mrk 509 (Kriss 
et al. 2000b; Kraemer et al. 2003; Yaqoob et al. 2003); NGC 3516 (Kraemer et al. 2002); 
NGC 3783 (Kraemer, Crenshaw, & Gabel 2001; Kaspi et al. 2001; Blustin et al. 2002; Gabel 
et al. 2003a,b); NGC 4051 (Collinge et al. 2001; Ogle et al. 2004); NGC 4151 (Crenshaw 
et al. 2000; Kraemer et al. 2001); NGC 5548 (Mathur et al. 1999; Crenshaw & Kraemer 
1999; Brotherton et al. 2002; Kaastra et al. 2000, 2002; Crenshaw et al. 2003; Steenbrugge 
et al. 2003). However, with the exception of Mrk 279 and some data taken for the long- 
term campaign on NGC 3783, these observations have been performed at different times 
in the UV and X-rays. The high degree of variability in AGN continua complicates self- 
consistent photoionization modeling of the intrinsic absorption from these non-simultaneous 
data, preventing firm conclusions about the nature of the relationship between UV and X-ray 
absorbers. 

The Seyfert 1 galaxy NGC 7469 (z = 0.01639, de Vaucouleurs et al. 1991) has been 
studied extensively in the optical, ultraviolet, and X-ray regimes (Wanders et al. 1997; Collier 
et al. 1998), and has been the subject of simultaneous UV and X-ray variability studies 
(Kriss et al. 2000a; Nandra et al. 1998, 2000). The intrinsic absorption in NGC 7469 has 
been studied previously in the X-ray with XMM-Newton (Blustin et al. 2003) and in the 
UV with FUSE (Kriss et al. 2003). However those observations were separated by one year, 
making the campaign presented here the first set of simultaneous, high-resolution UV and 
X-ray spectral observations. These authors found two primary components in the UV, with 
outflow velocities of —569 and —1898 km s -1 . The low velocity component is identified with 
the high-ionization phase of the X-ray absorption, also responsible for the X-ray emission. 
Its partial covering fraction suggests its position coincides with that of the BELR. The high 
velocity UV component has no associated X-ray absorber. Likewise, the low-ionization phase 
of the X-ray absorption appears to have no UV counterpart. 

Here, we present the simultaneous UV and X-ray observations and describe the proper- 
ties of the AGN outflow. We present photoionization models of the absorbing gas and draw 
conclusions about its geometry. 
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2. X-ray: Chandra 
2.1. Data 

We observed NGC 7469 on 2002 December 12-13 for 150 ksec with the High Energy 
Transmission Grating Spectrometer (HETGS) on the Chandra X-Ray Observatory. The 
HETGS 1 provides resolving power up to ~1000 from 0.5 keV to 8.0 keV through the use 
of two gratings: (1) the Medium Energy Grating (MEG) which gives good spectral coverage 
from 0.5 keV to 5.0 keV (2.5-25 A); and (2) the High Energy Grating (HEG), which is 
optimized for the 0.8-8.0 keV (1.6-15 A) spectral region. 

We extracted the spectra using the standard CIAO v3.0.1 pipeline. In order to maximize 
signal-to-noise, we co-added the ±1 order MEG and HEG spectra and binned to 0.02 A for 
analysis. Coarser 0.06 A bins were used to determine the continuum model. We used 
ISIS v.1.2.1 2 (Houck & Denicola 2000) for spectral analysis. See Table 1 for a summary of 
the Chandra observation details. 



2.2. Analysis: Continuum Fit and Spectral Features 

We fit the coarsely-binned MEG and HEG spectra in the 0.5-8.5 keV range with a power- 
law plus blackbody continuum modified by absorption from a Galactic hydrogen column of 
4.8 x 10 20 cm" 2 (Dickey & Lockman 1990). We used the xspec models tbabs (Wilms, Allen, 
& M c Cray 2000) and diskbb (see, eg. Mitsuda et al. 1984; Makishima et al. 1986), respectively 
to model the absorption and black-body components. The photon index of the power-law 
spectrum is T = 1.7910 ± 0.0009, with a normalization of A = (6.97 ± 0.10) x 10" 3 photons 
cm -2 s -1 keV -1 at 1 keV in the emitted frame, and the temperature of the blackbody is 
120±5 eV, with overall x 2 /dof = 1.12. As Blustin et al. (2003), we note that this continuum 
model is simply a convenient parametrization that allows us to fit the emission and absorption 
features in the Chandra spectrum, it is not meant to be a physical model of the continuum. 
The errors quoted are 90% confidence limits. The HETG spectrum, the continuum fit, and 
its components, including a Gaussian fit to the Fe Ka line, are shown in Figure 1. We also 
included additional multiplicative model components to account for calibration effects. We 
use: (1) an ad hoc correction for ACIS pileup effects which also reduces the effect of the 
iridium M-edge (Marshall et al. 2004a); (2) a correction to bring the frontside-illuminated 



1 For further details on Chandra instruments, see http://cxc.harvard.edu/proposer/POG. 
2 http://spacc.mit.cdu/ASC/ISIS/ 



- 5 - 



chip QEs into agreement with the backside-illuminated chip quantum efficiencies (Marshall 
et al. 2004a); and (3) an ACIS contamination correction (Marshall et al. 2004b). 

In regions of overlap, we require a spectral feature to be identified in both the HEG 
and MEG spectra in order to qualify as a bona fide absorption or emission line. Because the 
resolution of the HEG spectrum is 2 x that of the MEG spectrum, we use the HEG spectrum 
alone to measure the centroid equivalent width(flux) and FWHM of an absorption(emission) 
line in the 1.6-12.6 A region. At longer wavelengths, the S/N in the HEG spectrum falls 
below ~3 and we use the MEG spectrum. We also require reasonable oscillator strengths 
for quoted absorption features. For example, if we detect a Ly/3 feature from a particular 
species, we require its Lya transition to be detected as well. The HEG and MEG spectra 
are shown in detail, with the absorption and emission lines we identify and describe below, 
in Figures 2-4. The continuum fit shown in these figures is the global fit discussed above. 
Because of the low count rate in the ACIS pixels, we derive the flux errors per pixel from 
a N 1 + yj N + I where N is the count rate (Gehrels 1986), in order to approximate 
Gaussian statistics. For comparison we also show the XMM/RGS spectrum from Blustin et 
al. (2003), in the regions of overlap with the Chandra/BETG data, 6 — 25 A, in Figures 3 
and 4. 

We find several absorption features in the Chandra/EETG spectra, and we fit Gaussians 
to the unresolved profiles to derive the equivalent widths, the line centroids, and the velocity 
widths. We find Lya lines of the hydrogen-like ions Ne X and Mg XII, and a marginal 
Si XIV feature (Figure 5) and the Hea lines of Ne IX, Mg XI, and Si XIII (Figure 6). We 
also detect Ne IX He/3 with an outflow velocity and line width consistent with the Hect line. 
Of these features, Blustin et al. (2003) find only Ne IX He/5 in their XMM/RGS spectrum of 
NGC 7469. However we note that the resolution of the RGS is lower than that of both the 
HEG and the MEG, by factors of ~2.5 and ~5 respectively. Also, the Si XIV line lies at the 
high energy edge of the RGS spectral range. We also find three iron features, Fe XVII A15.0, 
Fe XXI A12.3, and Fe XX A12.8. No Fe XVII absorption is not found in the XMM/RGS 
data, but there is an ambiguous feature at the position of Fe XXI. The Fe XX feature falls 
on a gap in the data. Conversely, Blustin et al. (2003) report several absorption features 
that we do not find, or find only marginal evidence for, in the Chandra spectra: O VIII Lya 
and Ly (3 which we discuss further below, O VII 3-1 (r) (18.6 A), and Ne VII Beet (13.8 A). 
The other lines reported by these authors lie in regions where the HETGS effective area is 
low, or outside its spectral range altogether. 

The parameters of the absorption line fits are summarized in Table 2, as well as those for 
the emission features discussed below. Along with the equivalent widths, line centroids, and 
line widths, we tabulate the statistical significance of each feature defined by the detection 
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threshold. This is defined as the equivalent width divided by the the la equivalent width 
error, calculated as a function of wavelength using the spectrum flux and error boxcar- 
smoothed over ~ 2.5 x the FWHM of the line spread function (Bechtold et al. 2002). We list 
all the absorption and emission lines we detect at greater than ~ 2.5a in Table 2, including 
those for which we could make no identification. The low confidence level we set leaves open 
the possibility that some of the features in Table 2 are spurious. 

We find a marginal narrow Lja emission feature redward of the Ne X absorption, shown 
in the top panel of Figure 5. The significance of both the Ne X absorption and emission lines 
in Table 2 is low, 1.3a, and they are not detected in the XMM/KGS spectrum of Blustin et 
al. (2003). However, we do detect them at 4.1a and 3.6a in the MEG spectrum. 

We see some evidence of O VIII absorption blueward of, and possibly filled in by the 
Lja emission line, as reported by Blustin et al. (2003), but we do not find O VIII Lj/3 or 
Ly7 absorption in the Chandra spectrum (Figure 7) as Blustin et al. (2003) did. We do not 
find O VIII Ly/3, Ly7, or Lj5 emission corresponding to the prominent Lja emission line, 
but we do find significant emission features near the expected positions of the Ly-5, Ly-6, 
and Ly-7 lines. However, because the implied velocities of these features are not entirely 
consistent with the Lja line, these are marked as highly tentative identifications in Table 2 
and in Figure 7. These lines are not seen in the XMM/RGS spectrum. 

From the S/N in the spectrum at the expected position of O VIII Lja absorption, we 
place a 3a upper limit of ~30 mA on its equivalent width, which implies iV(0 VIII) < 
8.8 x 10 17 cm -2 for b ~ 100 km s -1 . We find no O VII Hea absorption, leading to an 
upper limit on the equivalent width of 80 mA, or N(0 VII) < 1.4 x 10 19 cm -2 , also for 
b ~ 100 km s _1 . This choice of Doppler parameter derives from the identification of this 
unresolved X-ray absorption with Component 1 in the UV spectra, which we model with nine 
subcomponents with FWHM=61 km s _1 , as described in Section 3.2.3. If we instead assume 
b = 200 km s _1 , the O VIII and O VII curves of growth give JV(0 VIII) < 5.1 x 10 16 cm" 2 
and JV(0 VII) < 3.4 x 10 18 cm" 2 , while b = 500 km s _1 gives iV(0 VIII) < 2.8 x 10 16 cm" 2 
and N(0 VII) < 5.8 x 10 16 cm -2 . These dramatically lower O VII column densities do not 
change the conclusions we draw from the photoionization models in Section 4. 

An emission feature consistent with O VII forbidden at —200 km s _1 is shown in the 
top panel of Figure 8. At ~110 mA, it is stronger than reported by Blustin et al. (2003). 
The corresponding intercombination line, however, is highly uncertain if it is present at all, 
as shown in the bottom panel of the figure. Blustin et al. (2003) report the forbidden to 
intercombination ratio f/i ~ 2 for O VII in the XMM data from 2000. If that ratio held for 
this observation, the S/N in the Chandra spectrum is such that we would only expect to find 
the intercombination line at ~ 2.5er confidence. We find a significant emission feature near 
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the expected position of the O VII radiative recombination continuum, listed in Table 2. 
The emission feature in the top panel of Figure 8 consistent with O VII(f) at —980 km s" 1 
is not statistically significant. Finally, we identify a significant emission feature at 13.91 A 
observed with Ne IX(f) at —200 km s _1 , although, again, we find no corresponding resonance 
or intercombination lines. 

A strong Fe Ka emission line is present in the Chandra spectrum. BeppoSAX obser- 
vations (De Rosa et al. 2002) indicated the presence of a narrow emission component in 
addition to the broad emission from an accretion disk. However, Blustin et al. (2003) found 
that the Fe Ka emission line in their XMM spectrum was best fitted by a single narrow 
Gaussian. Using the continuum defined above and fitting only the 1.5-2.5 A region of the 
HEG spectrum, we find a single Gaussian with flux 3.9 ± 0.7 x 10~ 5 photons cm -2 s -1 , line 
width 6310 ± 1580 km s -1 (FWHM), and rest-frame energy 6.39 ± 0.01 keV to be the best fit 
to this line with x 2 /dof= 62.5/57. Adding a second, narrow Gaussian does not improve the 
fit. These parameters are consistent with the XMM results. The width is better constrained 
by these Chandra data, however, and its value is consistent with the UV emission lines we 
fit in the STIS spectra, particularly the intermediate velocity width lines of Lya, C IV, and 
N V. These are discussed in the next section. We note that Nandra et al. (2000) found a sig- 
nificant correlation of the flux in the Fe Ka line and the 2-10 keV continuum flux, indicating 
that at least some of the flux arises within 1 light-day of the continuum source. 



3. Ultraviolet: FUSE and STIS 
3.1. Data 

We observed NGC 7469 with FUSE using the 30" x 30" low-resolution aperture. See 
Table 1 for the details of the observations. We combined the two exposures taken on 2002 
December 13 for a total exposure time of 6948 sec. Because of visibility constraints imposed 
by the coordinated observation, there was insufficient time to align the SiCl, SiC2, and LiF2 
channels for the FUSE exposures. Our method for correcting for the worm feature in the 
LiFIB segment uses the shape and flux level of the target in the LiF2A channel, therefore 
we were not able to perform the correction to derive a reliable spectrum over 1093-1182 A. 

We used the standard FUSE calibration pipelines (CALFUSE v2.2.3, see Sahnow et al. 
2000), to extract the spectra and to perform background subtraction and wavelength and 
flux calibrations. We combined the spectra from each detector and binned the spectra to 
0.05 A bins, improving the S/N while preserving the full spectral resolution, ~20 km s 1 . 

To correct the wavelength scale of the FUSE spectra for zero-point offsets due to place- 
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ment of the target in the 30" aperture, we applied a 0.0638 A shift to bring the wavelength 
scale into agreement with the heliocentric spectrum observed by FUSE in 1999 described 
by Kriss et al. (2003). These authors established the zero point of the wavelength scale by 
comparing the positions of interstellar absorption lines such as Ar I, Fe II, O I, and H 2 to the 
velocity of the Galactic 21 cm emission. We estimate that the residual systematic errors in 
the FUSE wavelength scale are on the order of 15 km s _1 , and that the flux scale is accurate 
to -10%. 

Simultaneously with the FUSE observations, we obtained observations of NGC 7469 
with the STIS FUV MAMA over 5 orbits (13013 sec) through the 0.2" x 0.2" aperture using 
the medium-resolution E140M echelle grating covering 1150-1730 A. We also obtained 9 
orbits (22810 sec) of STIS data on NGC 7469 over two visits on 2004 June 21-22 using the 
same set-up as above. See Table 1 for other details of the STIS observations. We used 
CALSTIS v2.13b to process the spectra. As described above, we measure the positions of 
lines from low- ionization species, Si II, N I, Si III, C II, C II*, Si IV, Fe II, and Al II in 
the STIS spectrum, and compare these to the mean heliocentric velocity of the H I 21- 
cm emission. We find that no correction to the wavelength scale derived by the CALSTIS 
pipeline is necessary. Using a white dwarf spectrophotometric standard, we applied a small 
correction to the flux in the STIS spectra to correct for residual echelle ripple structure and 
changes in the MAMA detector sensitivity not accounted for in the current STIS pipeline. 
Gabel et al. (2005) will present a full description of this process. 

Finally, we discuss the origin of a depression in the STIS spectrum at 1418.8 A. The 
expected position of redshifted Si IV at rest with respect to the AGN is 1416.61 A, so this 
feature is inconsistent with Si IV absorption due to intrinsic, outflowing gas. Also, this 
feature is shallow and broad, and qualitatively unlike other intrinsic features. We attribute 
this feature to an instrumental effect, namely a shadow on the MAMA detector. The STIS 
FUV MAMA employs a field electrode, or repeller wire, to direct any electrons emitted in 
directions away from the microchannel plate detectors back into the channels. The repeller 
wire shadow runs from (x, y) = (1,994) to (1048,1134). The feature appears in order 104, 
which is centered at y — 1030.7 in pixel space, and runs from y = 1038 to y = 1014, 
covering wavelengths 1414.3 to 1429.9 A. The wire crosses the dispersion axis at an obtuse 
angle, thus the feature is broad in the extracted spectrum. Because the shadow appears in 
slightly different locations depending which optical element is present in the spectrograph, 
it is not divided out in the flat field, but simply flagged in the data quality array. A broad 
depression in the flux in STIS/E140M spectra has been reported by others (Kraemer et al. 
2001; Leitherer et al. 2001), and a very similar feature is visible in an E140M spectrum of 
Mrk 279 (Gabel et al. 2005). Therefore, we conclude that this is not a feature intrinsic to 
NGC 7469. 
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3.2. Analysis 

3.2.1. Continuum and Emission Lines 

We used the IRAF 3 task specfit (Kriss 1994) and spectral regions unaffected by ab- 
sorption features to fit the continuum and emission lines in the FUSE and STIS spectra of 
NGC 7469. This routine is an interactive spectral fitting tool that uses x 2 minimization tech- 
niques to derive the best fitting parameters for a user-specified model. For the continuum, 
we fit a power law of the form j\ = /iooo^~ a to the 2002 FUSE and STIS spectra simulta- 
neously, where the normalization, /iooo> is the flux at 1000 A, and we include the Galactic 
extinction law of Cardelli, Clayton & Mathis (1989) with R v = 3.1 and E(B-V)=0.069 
(Schlegel, Finkbeiner, & Davis 1998) in the continuum fits. The best fit parameters are 
a = 1.082 ± 0.021 and / 100 o = 8.18 ± 0.06 x 10" 14 ergs s^ 1 cm" 2 A" 1 . For the 2004 STIS 
spectrum, we fit fx = 3.25 ± 0.01 x i(r 14 A- L086±a007 ergs s" 1 cm" 2 A -1 . We fit Gaussian 
profiles to the emission lines in all the spectra, the most prominent of which are broad and 
narrow components of O VI in the FUSE band, and Lya and C IV in the STIS band. The 
Lya profile requires four Gaussian components for a reasonable fit- broad, intermediate and 
two narrow components. The C IV emission line requires one broad, one intermediate, and 
one narrow component for a good fit, and we fit intermediate width and narrow components 
to the N V emission. These emission lines are plotted in Figures 9 and 10, and the fit 
parameters are tabulated in Table 3. 



3.2.2. Interstellar and Intergalactic Absorption Features 

There are 4 velocity components of 21-cm emission along the NGC 7469 sightline, includ- 
ing one high velocity component associated with the Magellanic Stream at —333 km s _1 (Wakker 
et al. 2001). The strongest component of the interstellar absorption, with strong molecular 
hydrogen absorption (Wakker et al. 2003), lies at zero velocity. We divide the ISM absorp- 
tion from the FUSE spectrum using a two-temperature model of the zero velocity absorption 
constructed from the FUSE spectral simulation code FSIM. Because the STIS spectra are 
much less crowded, we can use an empirical approach, modeling the ISM lines in the vicinity 
of intrinsic lines by fitting unblended lines of the same species and scaling by the appropriate 
oscillator strength ratios from published atomic data. 



3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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We identify an intergalactic absorption feature at z — 0.0134 from its strong Lya line 
in the STIS spectra. The qualitative shape of this feature is not as broad and shallow as the 
feature at 1419 A caused by the repeller wire shadow discussed above, and no such feature 
is expected at this position, or flagged near this position in the data quality array. 

Curiously, this intergalactic Lya lies between the two components of intrinsic absorption 
we discuss below. We consider this system intergalactic for the following reasons: (1) We 
find no O VI, N V, or C IV absorption corresponding to this system; and (2) The absorption 
does not vary between the 2002 and 2004 observations as the UV flux decreases by a factor 
of 2.5. This could because its absorption is highly saturated, but we rule out this possibility, 
even for a small covering fraction because no corresponding Ly(3 absorption is seen in the 
FUSE spectrum. Fixing the covering fraction to unity, we use specfit to derive N(R I) = 
6.7 ± 0.3 x 10 12 cm -2 and b = 52 ± 2 km s _1 for this intervening absorber. 

3.2.3. Intrinsic Absorption Features 

Kriss et al. (2003) reported two primary components of UV absorption, labeled Com- 
ponents 1 and 2, each with several velocity subcomponents. We use the IRAF task specfit 
to fit the absorption profiles as Gaussians in optical depth, specified by four parameters: (1) 
the optical depth at line center, (2) the line centroid, (3) the line width, and (4) the covering 
fraction. For doublets and members of a line series, we fix the central optical depths to 
the ratios required by the products of the wavelengths and the oscillator strengths. From 
these fits, we derive column densities and covering fractions for all subcomponents of the 
absorption complexes, Components 1 and 2. The errors on these parameters are determined 
by specfit by evaluating and inverting the curvature matrix around the final value of x 2 - F° r 
the fits to the absorption profiles, we fix the continuum and emission line fits to those derived 
from unabsorbed regions, as described above. 

Analysis of high S/N spectra by Arav et al. (2002), Gabel et al. (2003, 2005), and Scott 
et al. (2004) shows that the profiles of intrinsic absorption lines in AGN are best described 
using a velocity-dependent covering fraction. The FUSE and STIS data for NGC 7469, 
however, are too noisy for a reliable application of this technique. In addition, because we 
lack the full Lyman series from the FUSE data and because the C IV and N V doublets 
are weak in Component 1, as we discuss further below, we do not attempt such a detailed 
analysis of these data. The next-best approach, as used by Kriss et al. (2003), decomposes 
the line profile with multiple blended Gaussians. Properly fixed ratios of optical depth in line 
doublets enables us to constrain the covering fraction in the individual components. Kriss 
et al. (2003) found that the quality of the FUSE data were such that a single value of the 
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covering fraction adequately described each subcomponent. For consistency, and because 
our data have even lower S/N, we adopt the same approach here. 

In Component 2, we fit four subcomponents as did Kriss et al. (2003). However, in 
Component 1, we find blue and red wings of absorption in Lja, leading to a fit consisting of 
nine subcomponents. We found that our fits did not require any change in the covering frac- 
tion of the absorption in either component in either the 2002 or the 2004 epoch observations. 
We therefore fix the covering fractions to the 1999 values, 0.53 for Component 1 and 0.93 for 
Component 2. We fix the width of all subcomponents at 61 ± 4 km s~ x and 48 ± 5 km s _1 
(FWHM) for Components 1 and 2, respectively. The resulting outflow velocities and column 
densities in all the absorption subcomponents are listed in Table 4. For the purposes of our 
analysis, we use the summed column densities of all the subcomponents in Components 1 
and 2. The N(R I)-weighted outflow velocity of Component 1 is —562 ±11 km s" 1 and that 
of Component 2 is —1901 ± 3 km s -1 . These are plotted in Figures 11, 12, and 13 for all the 
species identified in intrinsic absorption in the FUSE and STIS spectra. 

The H I absorption is identified solely by the presence of the strong Lya lines in the 
STIS spectrum, seen in the top panels of Figures 12 and 13. In the 2002 FUSE spectrum no 
significant Ly/5 absorption is detected. For Component 1, this is consistent with the low H I 
covering fraction of the system. The Ly/3 in Component 2 lies at the position of a strong 
interstellar C II line. See Figure 11. 

The O VI doublet is strong in both components. As noted by Kriss et al. (2003), 
the strengths of the doublet lines in Component 1 are equal, or nearly so as best as can 
be determined by the signal-to-noise ratio in the 2002 spectrum, which results from an 
exposure ~5.5 times shorter than the 1999 spectrum. Figure 14 shows both lines of the 
O VI doublet for Components 1 and 2. It is clear from this figure that in Component 1, 
the flux ratio in the two line profiles is consistent with one within the uncertainties. The 
same is true, however, for Component 2 in this spectrum. We performed the initial fit the 
O VI absorption in Component 1 by fixing the velocity offsets and column density ratios of 
all the subcomponents with respect to the strongest subcomponent, Id, to those measured 
for the strong Lya absorption in the STIS spectrum. We found that allowing the velocity 
offsets of the subcomponents to vary from the fiducial Ly« values did reduce the \ 2 of the fit 
enough to justify the extra degrees of freedom, and the values in Table 4 reflect the improved 
fit. Though we are able to measure the O VI column density in Component 1 formally, by 
summing all the subcomponents, we treat the measurement as a lower limit, since the higher 
S/N 1999 FUSE spectrum showed that this component is highly saturated. This limit is 
shown in Figure 17, which we discuss further below in Section 3.2.4. We note the possibility 
that the O VI column density is different in these 2002 observations than in the 2000 epoch, 
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because the UV continuum level has decreased and because iV(H I) increased from its 1999 
value. We see no evidence for damping wings in the O VI profile, so an increase in N(0 VI) 
from 1999 would imply that the line is only further saturated and our ability to measure the 
column density precisely has not improved. We do not assume saturation for Component 2, 
because in several velocity bins, the flux ratio between the blue and red profiles is significantly 
greater than one, and again, because the 1999 spectrum indicates that the absorption is only 
nearly, not fully saturated. 

We identify prominent N V and C IV doublets associated with Component 2 in the STIS 
spectra. These species, including Lja absorption, are all identified in the 1996 epoch data 
from the HST/F&int Object Spectrograph and from the International Ultraviolet Explorer 
presented by Kriss et al. (2000). However, we also identify weak but significant N V and 
C IV doublets associated with Component 1. These are shown in Figures 12 and 13. While 
the red lines of each doublet are not obviously present in the spectrum, the blue members 
of the doubles are visible in at least some subcomponents. We established the significance 
of these doublets by comparing the x 2 °f a fit with no N V and C IV to that of a fit which 
did include these doublets. We fit these systems by fixing the covering fractions and line 
widths as discussed above. We allowed the column density and line center of subcomponent 
Id to vary and fixed the velocity offsets and relative column densities of all the other other 
subcomponents to those of the O VI subcomponents. For N V, the improvement in the x 2 is 
dramatic in both the 2002 and 2004 STIS spectra: A% 2 = 152 and A% 2 = 343, respectively, 
for the extra 2 degrees of freedom. For C IV we find A% 2 = 22 and A% 2 = 201. From this 
we conclude that these absorption features are indeed present in the data. The results of 
all these absorption profile fits are listed in Table 4 and the summed column densities are 
plotted in Figure 17. 



3.2.4. Variability 

Using previous FOS and STIS data along with the data presented here, we plot the 
variation of the 1000 A flux and the UV spectral index in Figure 15. From the 1996 obser- 
vation to 2002, the 1000 A flux increased by 26% while the power law index became harder 
by a similar factor. Between 2000 and 2002, the 1000 A flux decreased by a factor of ~1.6, 
and by another factor of 2.5 between 2002 and 2004. During this time, the spectral index 
got marginally softer, as expected, but the change is not significant given the uncertainties. 

The intrinsic Lya and C IV A 1548 absorption in the 2004 STIS spectrum is shown in 
comparison with that in the 2002 spectrum in Figure 16. We plot the total H I, O VI, 
C IV, and N V column densities in Components 1 and 2 in Figure 17. The top panel of the 
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figure includes the upper limit on N(H I) from the nondetection of Ly/3 in the 1999 FUSE 
spectrum presented in Kriss et al. (2003), as well as the measurements from Lja in the 2002 
and 2004 STIS spectra. Also, in the second panel, we show the N(0 VI) measurements for 
Component 2 and the lower limits from the saturated absorption in Component 1 in the 
1999 and 2002 FUSE spectra. 

From Figures 16 and 17, we see that the absorption in Component 2 does not vary in 
most species, with the exception of an increase in iV(H I) between 1999 and the 2002/2004 
epochs. In Component 1, the H I absorption grows stronger, consistent with a change in 
the column density of a factor of ~ 4 from 1999 to 2002 and again by another factor of 2.3 
from 2002 to 2004, as the UV flux decreased by a total of a factor of 4. The lower limit on 
iV(0 VI) is consistent with no change since the 1999 observation. The C IV column density 
increased by a factor of 4 between the 2002 and 2004 observations, and the N V column 
density also increased by a smaller factor as the UV flux decreased by a factor of 2.5. Thus, 
both velocity components of the intrinsic absorption respond over time to changes in the 
ionizing flux. In Section 5, we discuss how we use this variability to place limits on the 
radial distances of the absorbers from the AGN. 



4. Photoionization Models 

We computed photoionization models for the absorbing gas similar to those of Krolik 
& Kriss (1995, 2001) using the XSTAR photoionization code (Kallman et al. 2000). In 
Figure 18, we show the source spectral energy distribution used in the models. For SED1, 
we used the UV and X-ray spectral slopes and normalizations fit to the 2002 epoch spectra 
of NGC 7469, described above in Sections 2.2 and 3.2.1. For comparison, we show SED2, 
the SED based on the 2000 epoch FUSE and STIS observations of NGC 7469, presented 
used by Kriss et al. (2003). As described by Kriss et al. (2003), we assume f u oc u^ 1 at 
long wavelengths for both SED1 and SED2. This spectrum breaks at 2500 A to f u oc z/~ - 92 
for SED1 and /„ oc u~ - 75 for SED2 to match the 2000 and 2002 epoch FUSE spectra, 
respectively. The X-ray data are matched by f u oc z/~ as for SED1 /„ oc z/~ a7 for SED2 
between 0.5 keV and 100 keV. 

We ran models in two two-dimensional grids of total absorber column density and ion- 
ization parameter, appropriate for Components 1 and 2: (log N, logU) = (19.5-20.75, —3.0- 
4.0) and (18.3-19.0, —5.3 — 0.3), respectively, assuming solar abundances (Grevesse, Noels, & 
Sauval 1996) in both cases. We then compared the column densities of the observed species: 
H I, O VI, N V, C IV in the UV data, and the H-like and He-like Ne, Mg, and Si observed 
in the Chandra data, as well as limits on O VII and O VIII to the values predicted by the 
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models in order to determine the properties of the two velocity components of the intrinsic 
absorption. 

For Component 1, we find that U — 1.0 is consistent with the measurements of N V 
and C IV, with the lower limit on the saturated O VI absorption, and with the upper limits 
on the O VII and O VIII absorption from their nondetection in the Chandra data, for our 
assumption of b = 100 km s _1 , but also for Doppler parameters as large as ~ 1000 km s -1 . 
For the total column density, we find logiV = 20.0 is a reasonable match to the absorbing 
columns. For Component 2, we find (log TV, U) = (18.6, 0.08) gives the best match with the 
absorption data. 

All of the absorption from H-like and He-like ions in the Chandra data is seen at velocities 
consistent with the outflow velocity of Component 1. However, the column densities we 
measure imply larger ionization parameters from the models, 8 < U < 12, implying log iV > 
20, indicating the presence of multiple ionization phases, as found in the warm absorber 
associated with NGC 3783 (Netzer et al. 2003; Krongold et al. 2003) from spectral fitting to 
photoionization models. The relative dearth of spectral features in our Chandra spectrum 
of NGC 7469 precludes such detailed modeling in this case. 

5. Discussion 

These new simultaneous UV and X-ray observations of NGC 7469 allow us to compare 
the column densities of various species observed in the absorbers with photoionization models 
in a fully self-consistent manner in order to infer the physical conditions in the outflow. From 
our multi-epoch UV observations, we use the variability of the ionizing continuum and the 
absorption to draw conclusions about the distances of the absorbing components from the 
central engine of the AGN. 

Kriss et al. (2003) reported U = 6.0 for Component 1. Although this value is consistent 
with our measurement of O VI in the 2002 FUSE data, and the limits on O VII and O VIII 
from the Chandra spectrum, it is incompatible with the detection of any N V and C IV 
absorption in the 2002 STIS spectrum. With the added leverage afforded by the column 
density measurements of these species, we now estimate U — 1.0 for this velocity component. 
This decrease in the ionization parameter is consistent with the decrease in the ionizing 
flux between the 1999 and 2002 observations, although U — 1 is somewhat smaller than 
expected for a constant total column density and a scaling of the number of ionizing photons, 
which decreases by a factor of ~2 for the SEDs assumed for the 1999 and 2002 epochs 
shown in Figure 18. The changes in the H I, C IV, and N V column densities between the 
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STIS observations in 2002 and 2004 are also consistent with photoionization changes in the 
absorber. We have no knowledge of the X-ray continuum at the 2004 epoch, but if we make 
the simple assumption that the X-ray flux decreases in direct proportion to the UV flux, the 
column densities of these species are in rough agreement with logiV = 20 and U = 0.4. We 
lack the constraining power of a measurement of N(0 VI) at this epoch. 

As Steenbrugge et al. (2003) found for the warm absorber in NGC 5548 and Ogle et 
al. (2004) found for the warm absorber in NGC 4051, we find a trend of increasing column 
density with increasing ionization parameter for the species we detect in Component 1 of the 
NGC 7469 absorber. In Figure 19 we show the ionization parameter at which the abundance 
of each species peaks versus the equivalent hydrogen column density for that species in 
Component 1 with a least-squares estimation of the power law fit. 

The total hydrogen column densities derived from the Fe XX and Fe XXI lines would 
imply that the absorber is Compton thick. These columns are derived, as all the points 
in Figure 19 are, from the measured line equivalent widths using the curve of growth and 
a Doppler parameter of 100 km s _1 . We use the fractional abundance of the ion from the 
photoionization models and the assumed solar abundance of each atomic species. There are 
two ways the points could move down on Figure 19: (1) a larger Doppler parameter, eg. 
b = 500 km s _1 would yield iV H (Fe XX) = 2 x 10 23 cm' 2 and JV H (Fe XXI) = 3 x 10 24 cm" 2 , 
however, we note that this would affect all the points in the figure; or (2) an enhanced iron 
abundance in the absorber, as seen in some AGN broad line regions (Hamann & Ferland 
1999) and accretion disks (Lee et al. 1999, Ballantyne & Fabian 2001, Bianchi & Matt 2002, 
Schurch et al. 2003). In any case, removing the Fe XX and Fe XXI points from the fit does 
not change the slope or the standard deviation of the regression significantly. 

To treat the lower limit on Nn(0 VI) and the upper limits on iVn(0 VII) and iVn(0 VIII) 
we used the IRAF emmethod algorithm appropriate for censored data to perform the linear 
regression. In order to determine if the regression solution is driven primarily by the low-U 
points in Figure 19, we performed the linear fit without the C IV, N V, and O VI points 
included. The standard deviation of the fit and of the individual fit parameters increases, 
but within this increased uncertainty, the changes in the fit parameters are not significant. 
This fit is shown by the dotted line in Figure 19. This relationship between U and Nn 
illustrated in Figure 19 indicates that, while U = 1.0 for Component 1 is broadly consistent 
with the UV data, the absorber is composed of gas with a range of ionization parameters 
and the absorbing column scales with ionization. 

The column density and ionization parameter we infer for Component 2 are similar to 
the values inferred by Kriss et al. (2003) for the 1996 FOS and IUE observations (Kriss et 
al. 2000a). This is consistent with a change in the ionization parameter due to a change in 
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the ionizing flux, as the number of ionizing photons from SED1 is comparable to the number 
inferred from the Kriss et al. (2000a) SED. (See Figures 15 and 18.) 

We find X-ray absorption in helium-like and hydrogen-like neon, magnesium, and silicon 
all at velocities consistent UV Component 1 to within 2a, although not with the same 
total column density and ionization parameter as the UV-absorbing species. This supports 
the conclusions of Kriss et al. (2003) and Blustin et al. (2003), that this component is 
associated with the X-ray warm absorber and is located closer to the central engine than 
Component 2. This is further corroborated by the small covering fraction and outflow 
velocity of Component 1, ~ 0.5 and —560 km s _1 , respectively, compared with the unity 
covering fraction and ~ —1900 km s _1 velocity of Component 2 with respect to the systemic 
redshift of NGC 7469. The small covering fraction of Component 1 may indicate full covering 
of the continuum source and partial covering of the BELR. This geometry is consistent with 
the larger size of the BELR with respect to the central, continuum emission source, and has 
also been inferred in the outflow in Mrk 279 by applying simultaneous fits to several Lyman 
series lines and to the C IV, N V, and O VI doublets (Scott et al. 2004; Gabel et al. 2005). 
The increase in the H I, N V, and C IV column densities in Component 1 with the decrease 
in the flux in the ionizing continuum also supports this interpretation. 

The only X-ray absorption feature we find in common with the results of Blustin et 
al. (2003) is the Ne IX He/3 line at a rest wavelength of 11.5 A identified in Component 
1. We report its equivalent width as 14 ± 5 mA, while Blustin et al. (2003) found 50 ± 20 
niA. The marginal decrease in the Ne IX absorption is consistent with a photoionization 
effect and the adopted SEDs shown in Figure 18 in the sense that the flux at the ionization 
potential of Ne VIII is 23% larger in the 2000 epoch SED (SED1) relative to the 2002 SED 
(SED2). The total column density in Component 1 implied by UV absorption features from 
the photoionization modeling discussed in Section 4, logiV = 20, is consistent with that 
found for the X-ray absorber by Blustin et al. (2003) from their XMM data. However, as 
noted above, and as illustrated by Figure 19, the hydrogen-like and helium-like absorption we 
find in the Chandra data imply a higher total column density indicating a range of ionization 
phases. 

We use the fact that we see variability in the absorption in Component 1, particularly in 
its H I absorption, and the time between the 2002 and 2004 observations to place a 540-day 
upper limit on the recombination time of the gas, where t rec = in e a rec ^i^-)~ x (Krolik & 
Kriss 1995). For T = 3 x 10 4 K, the hydrogen recombination rate is ~ 2.3 x 10~ 13 cm 3 s -1 . 
This gives a fairly weak lower limit on the density, n e > 45 cm -3 , assuming a neutral 
fraction of 5 x 10 -4 . Using the ionization parameter U = 1.0, this gives an upper limit on 
the radial distance of Component 1 of ~100 pc. Following the prescription of Blustin et al. 
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(2005) we estimate the maximum radial distance under the assumption that the absorption 
arises in a region of width less than or equal to the distance from the ionizing source: 
r < L ion C v (j)/{^NH)-, where C„(r) is the volume filling factor. For the parameters of the 
X-ray absorption derived from the 2000 epoch Chandra observation of NGC 7469 (Blustin et 
al. 2003), Blustin et al. (2005) derive C„(r) = 0.086% and r < 1.6 pc for the X-ray absorber, 
which we argue is the same as Component 1 of the UV absorption. Using the parameters 
derived from the UV absorption, and assuming the same volume filling factor as Blustin et 
al. (2005), we find a somewhat firmer upper limit on the distance of Component 1 from the 
central engine than we could place from the variability, r < 15 pc. If the UV absorption 
arises in higher density material embedded in the X-ray warm absorber (Krolik & Kriss 
1995, 2001), its volume filling factor may be even smaller than that estimated above, and 
the radial distance correspondingly smaller. 

Blustin et al. (2005) found from their comprehensive analysis of the warm absorber data 
in the literature that most absorbers in Seyfert galaxies are located at distances from the 
central engine consistent with the torus, indicating that the outflow may originate there in a 
scenario similar to that proposed by Krolik & Kriss (1995, 2001). In the case of NGC 7469, 
they found that the warm absorber analyzed by Blustin et al. (2003) fits into this picture. The 
minimum and maximum distances they estimate for the absorber indicate that it lies outside 
the radius of the BELR, 0.0042 pc (Wandel, Peterson, & Malkan 1999) and within the radius 
of the torus, ~ pc (Krolik & Kriss 2001) where L is the ionizing luminosity between 1 
and 1000 Ryd. Allowing for the lower ionizing luminosity relevant to the 2002 observations, 
a slightly larger column density weighted ionization parameter for the X-ray absorption, and 
a mean outflow velocity somewhat lower than measured previously by Blustin et al. (2003), 
U ~ 10 versus U ~ 5, and —580 km s _1 versus —800 km s _1 , we find a similar overall 
result. For both the UV and X-ray absorption in Component 1, we find r min >> tbelr, 
given r min = 2G Mbh / v out- We note, however, that because the formal uncertainties on the 
black hole mass from Wandel et al. (1999), M BH = 0.76±g;?| x 10 7 M , allow the black hole 
mass to be as low as 1.5 x 10 5 M , or less, making r min consistent with the absorber arising 
from an accretion disk wind. Blustin et al. (2005) propose the alternative explanation that 
the wind has not yet achieved the escape speed and may therefore be closer to the central 
engine than r min as estimated above. This way, the partial covering of the continuum and/or 
broad line source demanded in particular by the saturated O VI profile of Component 1 may 
be understood. 

The larger outflow velocity of Component 2 compared to Component 1 leads to a smaller 
lower limit on the distance of Component 2 from the central engine, ~ 0.02 pc. However, 
given the unity covering fraction of this component, it is unlikely that we can place it within 
the BELR, even given the uncertainty in Mbh noted above. We use the change in iV(H I) 
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in Component 2 between the 1999 and 2002 epochs to estimate its maximum distance from 
the ionizing source from setting the recombination time to be less than the time between 
these observations. For a neutral fraction of 5 x 10~ 4 , t rec < 1100 days implies n > 20 cm -3 . 
The best-fit ionization parameter for this component, U = 0.08, implies r < 600 pc. 

We conclude that there is still ambiguity in the radial distances of the absorbers from 
the central engine of NGC 7469, but that the general picture described in previous work 
(Blustin et al. 2003; Kriss et al. 2003) holds: Component 1 is a high density, high ionization 
absorber located near or within the BELR and Component 2 is an absorber of relatively 
lower density and ionization located in a region that is consistent with an origin in the AGN 
torus. 
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Table 1. Observations of NGC 7469 



Instrument 


ID 


Start Date 


UT 


Exp. (s) 


\J 1 lilt bill ti 




2002- 


•12- 


•12 


lO.O t .uo 




Chandra 


700586 


2002- 


■12- 


•13 


12:10:13 


69760 


FUSE 


C0900101 


2002- 


■12- 


■13 


07:07:39 


3596 


FUSE 


C0900102 


2002- 


•12- 


•13 


06:25:22 


3352 


STIS 


O6BN01010 


2002- 


•12- 


•13 


06:28:04 


2233 


STIS 


O6BN01020 


2002- 


■12- 


■13 


07:51:30 


2695 


STIS 


O6BN01030 


2002- 


•12- 


■13 


09:27:36 


2695 


STIS 


O6BN01040 


2002- 


■12- 


•13 


11:03:42 


2695 


STIS 


O6BN01050 


2002- 


•12- 


•13 


12:39:47 


2695 


STIS 


O8N501010 


2004- 


•06- 


•21 


17:55:10 


1940 


STIS 


O8N501020 


2004- 


•06- 


•21 


19:14:52 


2870 


STIS 


O8N501030 


2004- 


•06- 


•21 


20:50:51 


2290 


STIS 


O8N501040 


2004- 


•06- 


•21 


22:26:50 


2870 


STIS 


O8N501050 


2004- 


•06- 


-21 


00:02:50 


2870 


STIS 


O8N502010 


2004- 


•06- 


•22 


17:55:3 


1940 


STIS 


O8N502020 


2004- 


•06- 


-22 


19:14:49 


2290 


STIS 


O8N502030 


2004- 


•06- 


-22 


20:50:49 


2870 


STIS 


O8N502040 


2004- 


•06- 


•22 


22:26:49 


2870 
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Table 2. Spectral Features in Chandra Spectrum of NGC 7469 



^vac 


Feature 


Equiv. width Flux 1 


Signif. 3 


Av 3 


FWHM 




(A) 


(mA) 


a 


(km s _1 ) 


(km s _1 ) 


Absorption 


6.18 


Si xiv Lya 


5±2 


2.3 


-540 ± 260 


620 ± 690 


6.65 


Si xill Hea 


11 ±6 


5.0 


-600 ± 330 


1350 ± 900 


8.42 


Mg XII Lya 


10 ±3 


2.6 


-540 ± 100 


500 ± 160 


9.17 


Mg XI Hea 


11 ±5 


3.3 


-620 ± 130 


680 ±410 


11.55 


Ne IX He/5 


14 ±5 


4.4 


-910 ±180 


1020 ±520 


12.13 


Ne X Lya 


17±7 


1.3 4 


-530 ± 80 


430 ± 270 


12.28 


Fe xxi 


21 ± 11 


2.2 4 


-540 ± 240 


890 ± 570 


12.82 


Fe xx 


27 ±6 


5.0 


-290 ± 130 


1070 ± 340 


13.45 


Ne IX Hea 


24 ±9 


3.6 


-670 ± 180 


1050 ±490 


13.45 


Ne IX Hea 


25 ±6 


3.2 


860 ± 100 


760 ± 220 


15.01 


Fe xvil 


19 ±8 


2.7 


-690 ± 150 


650 ± 320 


Emission 




Fe Ka 


3.9 ±0.7 


17 


6.39±0.01 5 


6310 ± 1580 


2.7 




1.0 ±0.4 


6.5 





3410 ± 1600 


8.2 




0.8 ±0.4 


6.2 





1600 ± 940 


11.0 




1.1 ±0.5 


4.5 





590 ± 300 


12.13 


Ne x Lya 


0.7 ±0.5 


1.3 4 


430 ± 170 


480 ± 400 


13.2 




1.0 ±0.4 


4.5 





790 ± 480 


13.69 


Ne ix(f) 


1.6 ±0.7 


4.8 


-200 ± 290 


1250 ±690 


14.45 


VIII Ly-7? 6 


1.0 ±0.4 


4.5 


-180 ±170 


780 ± 410 


14.52 


VIII Ly-6? 6 


1.4 ±0.6 


5.9 


810 ±280 


1270 ± 660 


14.63 


VIII Ly-5? 6 


1.5 ±0.8 


6.1 


1470 ± 470 


1790 ± 1190 


15.4 




1.3 ±0.6 


4.4 





820 ± 450 


15.5 




1.6 ±0.6 


5.5 





790 ± 340 


16.77 


VII RRC 


2.0 ±0.9 


5.8 


-1300 ±230 


1020 ±530 


17.1 




2.3 ± 1.0 


6.2 





1080 ±550 


17.3 




1.5 ±0.8 


3.8 





720 ± 460 


17.4 




1.4 ±0.7 


3.1 





490 ± 260 


18.3 




1.9 ± 1.0 


3.4 





700 ±410 
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Table 2 — Continued 



Avac 


Feature 


Equiv. width Flux 1 


Signif. 3 


Av 3 


FWHM 




(A) 


(mA) 


a 


(km s -1 ) 


(km s -1 ) 


18.97 


Vlll Lya 


3.0 ±1.0 


5.1 


50 ±90 


370 ± 340 


20.2 




2.4±1.5 


3.5 





620 ± 440 


22.09 


vil(f) 


6.1 ±2.5 


3.9 


-200 ± 90 


460 ± 220 



: In units of 10~ 5 photons cm -2 s _1 

2 Significance of feature from smoothed equivalent width threshold. See § 2.2. 

3 Velocity with respect to systemic redshift of NGC 7469, z = 0.01639 

4 Ne X absorption and emission are detected in the MEG spectrum at 4.1<r and 
3.6(7, respectively, and Fe XXI absorption is detected in the MEG spectrum at 3.4cr. 

5 Rest-frame energy of profile center in keV 

6 ID highly uncertain, no Ly/3, Ly7, or Ly<5 detected. See Figure 7 and discussion 
in § 2.2. 
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Table 3. Emission Line Fits to FUSE and STIS Spectra of 

NGC 7469 



Line 


A vac 


Flux 1 


Velocity 2 


FWHM 




(A) 




(km s _1 ) 


(km s _1 ) 


2002 FUSE 


VI broad 


1031.93 


3.77 ±0.31 


-137 ±180 


5907 ± 380 


VI broad 


1037.62 


1.88 ±0.16 


-137 ±180 


5907 ± 380 


VI narrow 


1031.93 


2.35 ±0.31 


-162 ±54 


1061 ±95 


VI narrow 


1037.62 


1.17 ±0.16 


-162 ±54 


1061 ±95 


2002 STIS 


Lya broad 


1215.67 


17.59 ±0.31 


-67 ±89 


13170 ± 130 


Lya int. 


1215.67 


17.13 ±0.49 


-202 ± 63 


4651 ± 131 


Lya narrow 1 


1215.67 


4.89 ±0.25 


-319 ±22 


1047 ± 55 


Lya narrow2 


1215.67 


3.26 ±0.52 


408 ± 50 


1586 ± 203 


N V int. 


1240.15 


5.22 ±0.05 


-107 ±33 


4651 ± 131 


N V narrowl 


1240.15 


0.66 ±0.06 


-107 ±33 


1047 ± 55 


Si II 


1263.31 


0.71 ±0.02 


-4 ±93 


2723 ± 174 


Si II 


1531.18 


2.09 ±0.11 


-82 ±53 


2170 ± 88 


C IV broad 


1549.05 


21.10 ±0.20 


224 ±51 


9371 ± 78 


C IV int. 


1549.05 


15.57 ±0.06 


-105 ±19 


3504 ± 32 


C IV narrow 


1549.05 


3.36 ±0.19 


-42 ± 12 


1118 ±47 


Fe II 


1608.45 


0.33 ±0.09 


795 ± 53 


2386 ± 315 


He II broad 


1640.50 


6.86 ±0.45 


-326 ± 297 


11041 ±706 


He II narrow 


1640.50 


0.89 ±24.41 


-326 ± 297 


1403 ± 188 


2004 STIS 


Lya broad 


1215.67 


4.84 ±0.06 


-1315 ±231 


13541 ± 144 


Lya int. 


1215.67 


11.66 ±0.12 


-129 ±23 


4774 ± 54 


Lya narrowl 


1215.67 


3.77 ±0.11 


-232 ± 16 


1137 ± 19 


Lya narrow2 


1215.67 


1.44 ±0.12 


449±9 


1566 ± 120 


N V int. 


1240.15 


3.84 ±0.09 


-38 ±15 


4774 ± 54 


N V narrowl 


1240.15 


0.42 ±0.03 


-38 ±15 


1137 ± 19 


Si II 


1263.31 


0.60 ±0.04 


-95 ± 117 


3428 ± 276 


Si II 


1531.18 


0.95 ±0.04 


-82 ± 53 


2844 ± 182 


C IV broad 


1549.05 


12.13 ±0.27 


381 ± 77 


9069 ± 149 
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Table 3 — Continued 



Line 


A vac 


Flux 1 


Velocity 2 


FWHM 




(A) 




(km s _1 ) 


(km s -1 ) 


C IV int. 


1549.05 


8.25 ±0.13 


-7 ±32 


3657 ± 38 


C IV narrow 


1549.05 


3.97 ±0.12 


-58 ± 18 


1458 ± 17 


He II broad 


1640.50 


2.79 ±0.75 


-326 ± 297 


10561 ±327 


He II narrow 


1640.50 


0.38 ±0.25 


-326 ±297 


1479 ± 913 



*Flux in units of 10 13 ergs cm 2 s 1 

2 Velocity relative to systemic redshift, z = 0.01639 
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Table 4. Intrinsic Absorption in FUSE and STIS 
Spectra of NGC 7469 



Feature 


Comp 


A vac 


N- 

1 ' ion 


Av 1 




# 


(A) 


(10 12 cm" 2 ) 


(km s _1 ) 


2002 


VI 


2a 


1031 


93 


19.2 ± 2.0 


-2014 ± 18 


vi 


2b 


1031 


93 


169 ± 17 


-1971 ± 2 


vi 


2c 


1031 


93 


690 ± 71 


-1897 ± 2 


vi 


2d 


1031 


93 


68.7 ± 7.0 


-1838 ± 18 


vi 


la 


1031 


93 


160 ± 15 


-730 ± 13 


vi 


lb 


1031 


93 


194 ± 18 


-628 ± 15 

V 7 \-J 1 _L ".7 


vi 


lc 


1031 


93 


439 ±41 


-594 ± 9 


vi 


Id 


1031 


93 


753 ± 70 


-546 ± 9 


vi 


le 


1031 


93 


309 ± 29 


-490 ± 12 


vi 


If 


1031 


93 


97.9 ± 9.1 


-483 ± 22 


vi 


le 


1031 


93 


93.4 ± 8.7 


-426 ± 34 


vi 


lh 


1031 


93 


76.8 ± 7.1 


-397 ± 23 


vi 


li 


1031 


93 


107 ± 10 


-320 ± 13 




2a 


1215 


67 


5.5 ± 3.1 


-2005 ± 18 


Lya 


9 V. 
ZD 


1 01 ^ 


O ( 


Q Q _l_ Q 1 
O.O Zt O.l 


1 QRD 4- 

— ryou zt oy 


Lya 


2c 


1215 


67 


200 ± 58 


-1897 ±2 


Lya 


2d 


1215 


67 


19.9 ±4.3 


-1841 ±5 


Lya 


la 


1215 


67 


23.3 ±3.5 


-736 ± 4 


Lya 


lb 


1215 


67 


28.1 ± 12.4 


-669 ± 15 


Lya 


lc 


1215 


67 


63.6 ±31.4 


-627 ± 14 


Lya 


Id 


1215 


67 


109 ±26.6 


-567 ±6 


Lya 


le 


1215 


67 


44.8 ±9.8 


-504 ± 13 


Lya 


If 


1215 


67 


13.5 ±9.1 


-460 ± 11 


Lya 


lg 


1215 


67 


14.3 ± 10.4 


-422 ± 22 


Lya 


lh 


1215 


67 


15.5 ±1.9 


-367 ±14 


Lya 


li 


1215 


67 


11.1 ±4.3 


-316 ± 19 


N v 


2a 


1238 


82 


6.7 ±0.3 


-2014 ±18 


N v 


2b 


1238 


82 


2.7 ±3.8 


-1971 ±2 


N v 


2c 


1238 


82 


242 ±11 


-1897 ±2 



-29- 



Table 4 — Continued 



Feature 


Comp 


'"vac 


^ * ion 


A 1 

Av 




# 


/ X N 

(A) 


(10 cm z ) 


n — 1 \ 

(km s J 


N v 


2d 


1 238 82 


24.1 ± 1.1 


— 1838 + 18 


N v 


1a 

X CI 


1 238 82 


11.4 ± 1.7 


—730 + 13 

1 UU _l_ X "J 


N v 


lb 

J. \j 


1 238 82 


38 ± 17 


—618 + 15 

UlU —1— X ?J 


N v 


lc 


1 238 82 


18 4 + 186 

1U. j: i 1U.U 


— 5Q4 + Q 

•J - J t: _l_ <J 


N v 


Id 


1 238 82 


53 5 + 8 

'.J . (J _l_ O . \J 


—546 + 9 

'J t: V _l_ iJ 


N v 


le 


1 238 82 


22 + 3 3 


—490 ± 12 


N v 


If 


1238 82 


6.6 ± 1.0 


—483 ± 22 


N v 


1 g- 
L & 


1 238 82 


7 + 10 


—426 ± 34 


N v 


lh 


1238 82 


5 4 + 08 


—397 ± 23 

"J U 1 _l_ '—' ' J 


N v 

IN V 


_L 1 


1 238 82 


7 6 + 11 

( . VJ _1_ X . X 


_Qon + 1 3 

UZjU l x*j 


C iv 

v. 1 V 


2a 


1 548 20 


5 3 + 02 

cj . - J _l_ VJ . 


—2014 ± 18 


C iv 


2b 


1548 20 

X OTTO . Zj w 


2 + 01 

_ . VJ _l_ VJ . X 


— 1971 ± 2 


C iv 


2c 


1 548 20 


190 + 8 


— 1897 + 2 


C iv 


2d 


1548.20 


18.9 ± 0.8 


-1838 ± 18 


C iv 


la 


1548.20 


1.6 + 0.4 


-730 ± 13 


C iv 


lb 


1548.20 


5.5 ± 1.4 


-618 + 15 


C iv 


lc 


1548.20 


2.6 + 0.7 


-594 ± 9 


C iv 


Id 


1548.20 


7.6 + 1.9 


-546 ± 9 


C iv 


le 


1548.20 


3.1+0.8 


-490 ± 12 


C iv 


If 


1548.20 


0.9 + 0.2 


-483 ± 22 


C iv 


lg 


1548.20 


1.0 + 0.2 


-426 ± 34 


C iv 


lh 


1548.20 


0.8 + 0.2 


-397 + 23 


C iv 


li 


1548.20 


1.1+0.3 


-320 ± 13 


2004 


Lya 


2a 


1215.67 


6.4 + 3.0 


-2005 ± 18 


Lya 


2b 


1215.67 


5.1 + 3.1 


-1959 + 38 


Lya 


2c 


1215.67 


191 + 36 


-1897 + 2 


Lya 


2d 


1215.67 


36.3 + 5.4 


-1841 + 5 


Lya 


la 


1215.67 


58.1 + 6.1 


-736 ± 4 


Lya 


lb 


1215.67 


0.09 + 5 


-669 ± 15 
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Table 4 — Continued 



Feature 


Comp 


A vac 


N ion 


A 1 




# 


(A) 


(lO 1 ^ cm z ) 


(km s J 




X L. 


1 91 67 

XZj X tj. U 1 


QfiA _|_ on 


_«97 _i_ i a 

VJZj i _l_ 11 




-LU. 




ZjZjU zl ou 


_c«7 _i_ a 




1 P 


1 91 ^ 67 

XZj X tj. U 1 


oo 7 _)_ i n o 


_KC\A + 1 Q 

uUt: _L Id 




1f 

X 1 


1 91 ^ 67 

X Zj X tj . \J ( 


9 8 ± 1 4 

t7 . O _1_ X . t: 


—460 ± 1 1 


T ,"\7'*'~V 

ijyix 




1 91 ^ 67 

XZj X tj. U 1 


Q1 7 + 4(1 


_499 _i_ 09 

TzZjZj _l_ ZjZj 


T .A/TV 


X 11 


1 91 ^ 67 

X Zj X xj . \J 1 


14 9 + 39 

XL<_/ l O.Zj 


—367 ± 1 4 




2i 

ZjI 


1 91 ^ 67 

XZj X tj. U 1 


7 q _i_ 9 « 

1 . C/ -XZ Zj. u 


_qi a _i_ i n 

O 1U X C/ 


N v 

IN V 


2a 


1 938 89 

X ZjOCj . (JZj 


6 ± 4 


—201 4 ± 1 8 

ZjUIt: l X O 


N v 

In V 


Zj U 


1 938 89 

X ZjOO . OZj 


8 Q + ^ Q 


_i q7i _i_ 9 

1 J 1 1 _l_ Zj 


N v 

In V 


2c 

ZjU 


1 938 89 

X ZjOO . OZj 


91 6 4- 1 ^ 

Zj X \J _l_ X o 


_i oq7 _i_ 9 

IOj 1 -J-l Zj 


N v 

IN V 


2d 

Zj kj. 


1 938 89 

X ZjOCj . (JZj 


21 Fi + i 


— 1 838 + 1 8 

IOlKJ _i_ xo 


N v 

IN V 


1 a 

-L CI 


1 938 89 

X ZjOCj . (JZj 




_7Qfj _|_ 1 Q 
1 UU _l_ xo 


N v 

In V 


J. U 


1 9^8 89 

XZjOO. OZj 


«4 n _i_ 90 « 




N v 

In V 


1 c 

X U 


1 938 89 

IZ/tJO. OZj 


41 n _i_ 9« n 


_cq4 _i_ q 


N v 

In V 


X LI 


1 938 89 

IZ/tJO. OZj 


sqn + 141 




N v 

In V 


1 p 


1 938 89 

IZ/tJO. OZj 


o« c: _j_ c o 

OU.u _I_ <J .O 


_4qn _i_ i o 

t:C/W _l_ XZj 


N v 

In V 


1 f 


1 9^8 &9 

XZjOO. OZj 


110-1-17 

ll.U ZL X . ( 




N v 

IN V 




1 938 89 

X ZjOO . OZj 


11 fi + 1 8 

X X . \J I X . (J 


—426 ± 22 

t:Zj\J 1 ZjZj 


N v 

In V 


1 h 

X 11 


1 938 89 

X ZiuO . OZj 


Q 1 4- 1 4 

C/ . X _l_ X . t: 


_oqq _i_ i q 

Octet ZC XO 


N v 

IN V 


1 i 


1 938 89 


1 9 fi + 9 


_Q9n _i_ oq 


L IV 


2a 


1 r in OH 

1548. 2U 


4.8 ± O.o 


— 2U14 ± 18 


C iv 


2b 


1548.20 


1.8 ±0.1 


-1971 ±2 


C iv 


2c 


1548.20 


173 ±12 


-1897 ±2 


C iv 


2d 


1548.20 


17.2 ± 1.1 


-1838 ±18 


C iv 


la 


1548.20 


6.6 ±0.6 


-730 ± 13 


C iv 


lb 


1548.20 


22.2 ±1.9 


-618 ±15 


C iv 


lc 


1548.20 


10.8 ±0.9 


-594 ± 9 


C iv 


Id 


1548.20 


30.9 ±2.7 


-546 ± 9 


C iv 


le 


1548.20 


12.7 ± 1.1 


-490 ± 12 


C iv 


If 


1548.20 


3.8 ±0.3 


-483 ± 22 
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Table 4 — Continued 



Feature Comp 


^vac 


iVion 


Av 1 




# 


(A) 


(10 12 cm" 2 ) 


(km s -1 ) 


C IV 


lg 


1548.20 


4.0 ±0.3 


-426 ± 34 


C iv 


lh 


1548.20 


3.1 ±0.3 


-397 ± 23 


C iv 


li 


1548.20 


4.3 ±0.4 


-320 ± 13 


1 Velocity 


with 


respect 


to systemic 


redshift of 


NGC 7469, 


z = 0.01639 (de 


Vaucouleurs et al. 1991) 
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Fig. 1 — Chandra HETG spectrum of NGC 7469, binned to 0.02 A, with best fit model 
including Fe Ka emission (red). Power law + iridium M-edge at 5.9 A (blue) and blackbody 
(magenta) components are shown separately. 
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Fig. 2. — HEG spectrum of NGC 7469 with flux errors with absorption (emission) lines 
identified below(above) the spectrum by tick marks with labels where ID known or "U" for 
unidentified lines. The continuum fit in Figure 1 is shown by the solid line. 
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11 11.5 12 12.5 

Rest Wavelength (A) 

Fig. 3. — HEG spectrum of NGC 7469 with flux errors with absorption (emission) lines 
identified below(above) the spectrum by tick marks with labels where ID known or "U" for 
unidentified lines. The continuum fit in Figure 1 is shown by the solid line. Also shown 
is the XMM/RGS spectrum obtained in 2000 by Blustin et al. (2003), offset by arbitrary 
factors for clarity. 
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Fig. 4. — MEG spectrum of NGC 7469 with flux errors with absorption(emission) lines 
identified below(above) the spectrum by tick marks with labels where ID known or "U" for 
unidentified lines. The continuum fit in Figure 1 is shown by the solid line. Also shown 
is the XMM/ RGS spectrum obtained in 2000 by Blustin et al. (2003), offset by arbitrary 
factors for clarity. 
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Fig. 5.— Hydrogen-like lines in Chandra/RETGS spectrum of NGC 7469, all from the HEG 
channel. The continuum level is the global fit shown in Figure 1. 
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Fig. 6.— Helium-like lines in Chandra/ HETGS spectrum of NGC 7469, Ne IX at v = -590 
and 860 km s _1 from MEG, Mg XI and Si XIII from HEG channel. The continuum level is 
the global fit shown in Figure 1. 
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Fig. 7. — Top left panel: VIII Lya emission, detected at 3.5a, and possible absorption 
in the Chandra MEG spectrum of NGC 7469. Middle and bottom left panels and top right 
panel show absence of Ly/?, Ly7, and Ly<5 features. Middle and bottom right panels show 
emission features listed in Table 2 tentatively identified as O VIII Ly-5, Ly-6, and Ly-7. The 
continuum level in all panels is the global fit shown in Figure 1. 
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Fig. 8. — VII forbidden emission and the expected position of corresponding intercombi- 
nation emission in the Chandra MEG spectrum. The continuum level is the global fit shown 
in Figure 1. 




Fig. 9. — The 2002 Dec FUSE spectrum of NGC 7469 (black) with errors shown in blue, 
and with continuum, emission line, and ISM absorption fits (red). 
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Fig. 10.— The 2002 Dec and 2004 Jun STIS spectra of NGC 7469 (black) with errors shown 
in blue, and with continuum and emission line fits (red). 
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Fig. 11.— Normalized 2002 FUSE spectrum shifted to rest frame of NGC 7469, with ISM 
fits (dashed lines) and Ly(3 and O VI intrinsic absorption. Vertical dot-dashed lines represent 
the column density-weighted mean velocity of all H I subcomponents in Components 1 and 
2, —562 ± 11 km s _1 and —1901 ± 3 km s -1 , respectively. 
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Fig. 12.— Normalized 2002 STIS spectrum shifted to rest frame of NGC 7469, with ISM fits 
(dashed lines) and Lya, N V and C IV intrinsic absorption. Vertical dot-dashed lines repre- 
sent the column density-weighted mean velocity of all H I subcomponents in Components 1 
and 2, —562 ± 11 km s" 1 and —1901 ± 3 km s" 1 , respectively. 
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Fig. 13.— Normalized 2004 STIS spectrum shifted to rest frame of NGC 7469, with ISM fits 
(dashed lines) and Lya, N V and C IV intrinsic absorption. Vertical dot-dashed lines repre- 
sent the column density-weighted mean velocity of all H I subcomponents in Components 1 
and 2, —583 ± 9 km s _1 and —1902 ± 3 km s -1 , respectively. 
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Fig. 14. — VI doublet in Components 1 and 2. The flux in the 1032 A line is shown by 
the solid bold lines, errors are shown by the dotted bold lines. The 1036 A lines are shown 
by the regular solid and dotted lines. 



-46- 



i 1 1 r 



15 



i 1 1 r 



1999 



n 1 1 r 



1996 



10 - 



2002 



o 5 



2004. 



o 
o 
o 



H 1 1 h 



H 1 1 h 



H 1 1 h 



1.3 



> 



1.2 



1.1 



i 



1 - 



o 



J I I L 



1 



J I L 



1000 2000 
MJD-50000 



3000 



Fig. 15. — Variation in the continuum flux and the UV spectral index of NGC 7469. 
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Fig. 16. — Comparison of absorption in Lja and C IV between 2002 and 2004. 
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Fig. 17.— Summary of H I, VI, C IV, and N V absorption intrinsic to NGC 7469 in 1999, 
2002, and 2004 epochs. 
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Fig. 18. — Spectral energy distributions used in photoionization models of the intrinsic 
absorbers in NGC 7469. SED1 is normalized such that it matches the 1000 A flux in the 
2002 epoch FUSE spectrum and the 2 keV flux in the Chandra spectrum. SED2 is normalized 
to the 2000 epoch FUSE and XMM-Newton data plotted, as described by Kriss et al. (2003). 
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Fig. 19. — Equivalent hydrogen column density versus ionization parameter at peak abun- 
dance of each ion with power law fit. The dotted line shows a fit with the C IV, N V, and 
O VI points omitted. 



